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Abstract 



r-| ! It is shown that amphtude-based, exact resummation tames the un-canceled IR diver- 

Q-i' gences at 0{al) in initial state radiation in QCD with massive quarks. Imphcations for 

Ph' precision predictions for LHC physics are discussed. 
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The era of precision QCD at the LHC, by which we mean 1% or better precision tags 
on the theoretical predictions, presents us with the extremely challenging task of proving 
that a given theoretical precision tag does in fact hold to that level. This means that 
all aspects of the standard formula for hadron-hadron scattering in perturbative QCD 
have to be examined for possible sources of uncertainty in the physical and technical 
precision components of any quoted total theoretical precision tag. In this connection, 
we note the standard practice of treating all quarks in the initial state as massless. It 
would be desirable to put an explicit error tag on this assumption by doing the respective 
calculations with the respective quark masses at their known [1] values and comparing 
the attendant predictions with their massless limits. This direct approach is however 
currently blocked by the pioneering results in Refs. [2,3], wherein it has been established 
that there is a lack of Bloch-Nordsieck cancellation at 0(0^) in the initial state radiation in 
massive QCD. Hence, even the b quark has to have zero mass in the initial state radiative 
corrections when one works to 0{a^) with n > 2. 

In what follows, we re-visit the results in Refs. [2] from the standpoint of recent 
progress [4, 5] in the resummation of large IR effects in the QCD perturbation theory, 
where we will focus on exact resummation methods with an eye toward rigorous control 
on any theoretical precision error budget that we may ultimately want to advocate. In 
this context, let us recall already the master formula that we have derived in Refs. [4]: 
using a. 2 ^ 2 + X hard process with multiple gluon (G) emission, q{pi) + q'{qi) — > 
q"{P2) +q"'iq2) + n{G) +X{px) in an obvious 4-momentum assignment notation, we have 
the differential cross section 

^^^^p ^ gSUMiR(QCD) — /" TT f ^^^^ f.iy(Pl+qi-P2-q2-Px-T. fc,)+goCD 

n=0 J j=l 

g , ^ d^p2 d^q2 d^Px 

* Pnl"^l; • • • ; "^nj n n n~ 
P2 (l2 Px 

where the hard gluon residuals Pn{ki, . . . , kn) and the infrared functions SUMir(QCD), -Dqcd 

are defined in Ref. [4] and we stress that the . . . , kn) are free of all infrared diver- 

gences to all orders in as{Q). (See especially Ref. [5] for explicit application of ([T]) to a 
real bremsstrahlung process.) Note that the hard gluon residuals Pn{ki, . . . , kn) have the 
structure [4] 

^n{kl, ...,kn) = [tri)M^^^^M^^'^]\lR.subtracted 

where the IR-subtraction is as given in Ref. [4] and f) is the initial-state color-spin density 
matrix so that the full quantum mechanical color effects are included in ([1]). For our 
initial state radiation (ISR) analysis, we take q and q' to be massive quarks of mass rUg, 
we take X to be a (QCD singlet) electroweak gauge boson to match the problem studied 
in Ref. [2] and we only compute ISR radiative effects in QCD. Let us then recall the 



^We do not employ here the resummation algebras that focus on the boundary of phase space [6,7], 
as these engender approximations which require arguments that are not opportune for our purposes. 
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pioneering result in Refs. [2]: working in the eikonal approximation and discussing the 
part of the cross section proportional to the color structure (here, H corresponds to the 
attendant hard sub-process) 

F, = C,{G)H:^{T,f-{T,),a (2) 

for the process — > V^*^ + X' where V^*"^ is our (off-shell) electroweak gauge boson 
and a, a are the colors of the quarks, the authors in Refs. [2] find the IR divergent result 



(PQ (rf- 4)327r2 V (3 J \(3 '1-/3 

, where g is the QCD coupling constant, H is the attendant hard sub-process factor dressed 
as in the color structure Fi, f] is the velocity of one quark in the rest frame of the other, 
and d is the dimension of space-time with c? > 4 to regulate the uncanceled IR divergence. 

is the invariant mass of the V*. This divergence is clearly non-Abelian in character as 
it vanishes for C2{G) = 0, where we define the gluon and quark representations' quadratic 
Casimir invariants respectively as usual: 

fijkfiji = C2{G)Ski 

{TiTijab = C2{F)Iab = Cplab, 

where fijk are the group structure constants. The result shows a clear lack of Bloch- 
Nordsieck cancellation at 0{al) and the standard approach is to set rriq = so that this 
uncanceled IR divergence vanishes as /3 — > 1 as one can see from ([3]). 

We point-out that the authors in Ref. [3] have analyzed the problem studied in Refs. [2] 
from a coherent-state Hamiltonian approach and have corroborated the result ([3]) with 
the added understanding that, in the coherent-state approach, the single pole divergence 
is converted into an unfactorizable, unspecified dependence of the respective coUinear 
singularities on the scale separating the attendant observable and un-observable gluon 
degrees of freedom inherent therein. This is again unacceptable and forces the use of 
rrig = for initial ISR for calculations at (9(a"), n >2. 

Here we propose an alternative approach. We look into the systematics of the analysis 
of the first paper in Refs. [2] . We see that one can represent the RHS of as the left-over 
real IR divergence which is uncanceled by the virtual IR divergence. In the language of 
the diagrams analyzed by Mueller's theorem [8] in the aforementioned paper, the RHS 
of ([3]) can be identified with a fraction Fnbn of the contribution of the real emission from 
the contribution of the diagrams equivalent to the diagram contribution (q-o) in Fig. 6 in 
the first paper in Ref. [2], which we reproduce here for definiteness in Fig. 1. To see this, 
let us recall the result of this last paper for the (q-o) diagrams (see Fig. 1) contribution 
to the differential cross section, removing the kinematic (note in this language the hard 
scattering factor is kinematic to the soft interactions under study here) and color factors: 
from the 5th equation on page 11 of the paper, we have the result, from the equivalence 
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of diagrams c and q and the equivalence of diagrams o and f in Fig. 6 of the paper (Fig. 1 
here) , 

'~°~]3^ J {k, + k', + le) {(3^kl - k^) {(3^kl - k'2 + le) {kl + e^Y ^' 

where we denote the 3-momentum by boldface letters, so that /c = k, and where the 
eikonal limit has been used in ([5]) as it was in respective paper in Ref. [2]. Because 
of this approximation, there is a spurious UV divergence in ([5]), which does not affect 
the IR regime. The authors in Ref. [2] therefore regulate this UV divergence with the 
factor e~^^/^^ for each would-be 3-space integral and then use dimensional methods [10] 
to isolate the IR divergence of interest; one obtains in this way the UV regulated result 
from ([5]) 

where the UV cut-off A is large compared to the soft scales in the problem, and here we 
have d = n + 1 io make contact with Q. When one adds the remaining contributions 
associated to the remaining graphs in Fig. 1, one sees from comparing ([3]) and ([6]) that 
the double pole term in ([6]) is canceled and that the fraction 



nbn 



(1-/5) (In (i±f)-2/3) 

— (7) 



In ('-±^\ 



of the single pole term is left over as the uncanceled IR divergence. 

The classic Landau-Bjorken-Cutkosky (LBC) analysis then allows us to determine the 
relationship between the real emission in the (q-o) diagram contribution and the single pole 
term on the RHS of ([6]). Specifically, upon doing the integral on the RHS of ([5]) over k' 
there are are two poles in the respective complex plane below the real axis, one at — k^ — ^6 
and one at —\J (S'^k'l — k' ±^ + ie, where here the energy of the /c'-gluon is just —Pk^ by the 
LBC rules in this eikonal exercise. The contribution of the former pole does not result 
in on-shell k' gluons . The LBC rules tell us that the regime 9^ = {0 < k'±^ < (3'^k'^} 
represents the regime wherein the fc'-gluon is actually on-shell here. Focusing on this 
regime, we see that we have the contribution 

1 r^su -2.. 1 



-(-2) ^//3^kl - k\' kz - .Jp^kl - k\' + ze 
2/cp- 



/32A;2-k2 kl + e^' 



where we have written the 2-space integration measure as 7rd{k'±^) by doing the respective 
angular integral. The integration over the latter measure can then be re-written, using 
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the fact that we only need the real part, 



' ' J Jo ^ ' Jem - k' , ik. + lef - W^kl - k' 



^j32kl - {k, + ze)2 - {/3^kl - k\^) 
1 2k. 



I3^kl - k2 kl + e2 



' ^(i-'kl - k\' {K + - {/3^kl - k'^') 
1 2k, 



3?—^ / d^k I Tld{k\^) 



-Til 1 
' J 


d^k j 












I d^k{- 





(3^kl -k^kl + e2 



^/P^kl - k'j_' 



2k, 



2K 



(9) 



where we again emphasize that the on-shell regime actually has /cq = —^kz < so that 
the real radiative contribution, by the standard LBC methods, has k^ > 0. If we integrate 
over the region kz > ^/e, it is clear that the RHS of the last equation has no real part 
as e — > 0. Thus, the real emission part of must arise from the regime < /c^ < ^/e. 
We treat the branch cuts for the logs by joining them between kzi = —ie/{l — (3) and 
kz2 = — + (3) and then we close the contour below the real axis as shown in Fig. 2 
to get the result, by Cauchy's theorem, 

1 2k 

dkz (- ln(A;, + ie- (3kz) + ln(A;, + ze + (3kz)) ^2^2 _ ^2 ^2 = 0' (1°) 

where we use the intrinsic freedom in the Feynman ie-prescription to take each such 
infinitesimal parameter independently to from above and the curve C is given in Fig. 2. 
We take here k± > If we denote the integrals over the i — th part of C by /j, i = 
1, ■ ■ ■ ,7, where the labels for these parts are defined in Fig. 2, then one can readily see 



^We use standard Lebesgue integration theory to conclude that the order of integration does not 
matter so that fixing k± and integrating over kz first, which means that the Hmit e —^ will always 
give us k± > ^/e, followed by integration over the full range of fc_L, when all integrals are finite by 
regularization where necessary, can not affect the final result. Alternatively, the reader can check that 
with the regularization we use, if one does the attendant integral over < fcj^ < ^/e, the respective result 
will vanish for e ^ 0. 
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that (for example we may set e = e 2 ) we have 

/v^ 1 2k 

7 



1=2 

We now treat the integrals /j, i = 2, ■ ■ ■ , 7 in turn. 

For J2, use the change of variable kz = -y/ee*^, for > > — |. Then, we get 

1 2kz 
(3^kl -\i^kl + e2 







— ivrln 



l + [3 



1-/9; (-k^^) 



2\ ■ 



1 ' idOkz (- ln(A;^ + ie - /3k^) + ln(fc^ + ie + (3k^)) 
Jo 

2 / (- ln(l + ln(l + /?)) — ^ (12) 

— k I 



For I3 it is enough to use the change of variable k^ = —iy to see that it is pure real 
so that it will not contribute to the imaginary part of Ji via (ITOl) and only this part of Ii 
is needed in extracting the real emission part of the RHS of (Q. 

For I4 we see from passing around the lower branch point in Fig. 2 that the respective 
imaginary contribution is just 

^ .... 1 2k, 



/1+/3 
dkz (— Tri) 



27ri In 



l + 1 



For Is, we see by the change of variable kz = —iy that it is pure real and does not 
contribute to the imaginary part of Ji via (fTUjl . 

For Iq, we get the result 

Ie = 7iiRes{-ie) 

since e/e ^ when e ^ 0. 

Finally, for J7 the change of variable kz = —iy shows that it too is pure real and does 
not contribute to the imaginary part of Ii via (fTOj) . 
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The net result is that we arrive at 



i^h = -{2TTi - vrij^ln (j^f ) 



In 



When we introduce the RHS of f|T5|) into iQ we get the result 

A_.U,.a,rad. = ^ (^In ([^)) (16) 

where we explicitly indicate that this is the real emission contribution by the subscript 
real rad.. Using the UV regulator employed in the first paper in Refs. [2], we see that the 
integral over k_|_ in (fT6l) can be written as 

(17) 

We regulate the infrared divergence by analytic continuation to n dimensions to get 



2^UV reg. — / ,2 



UV reg.jIR reg. 



k2 

k2/A2-pk2 



dp I d''k6{k,)e 



2 n n-3 

= ^(A') 2 . 

n — 6 

Introducing this last result into f|T6|) . we get 

A-ok.eal rad., UV reg. - J, V 1^ J J ' ^ ^ 

which shows that the real emission part of Aq^o saturates its single IR pole contribution. 

Isolating the divergent single pole IR term in f|T9l) we may now re-write the pioneering 
result of Refs. [2] as follows: the uncanceled IR singular contribution to the respective 
differential cross section is 

c/3g 647r 



da —g'^H r-, . , 

jS/n ~ — fi~-^n6n^iJ-o|9^,real rad., IR pole part; (20) 



where from (fT9l) we have 



_ 47r^ / 1 /I . 

^(j-o|sH,real rad., IR pole part — \2{n — 3) 1 ]^ _ ^ ) ) ' ( ^ 

6 



We note that the resuh in ( |T9l) agrees with the single pole term in ([6]) and with ( pTl) up 
to finite terms. 



From the result (l20l) we can now see how the theory of exact, amplitude resummation 
may impact the conclusions of Refs. [2] . We apply the formula in ([1]) to the real emission 
process in A^.^l^^, following for example the steps given in in Ref. [5]. We stress that 
we apply the resummation only to that fraction, Fnbn, of the real emission that has the 
uncanceled IR singularity in ( l20i) . The remaining 1 — Fnbn is not resummed because it 
is canceled by the sum of the remaining contributions associated with the diagrams in 
Fig. 1. We get, in this way, the result 

-^nfen^q— o I iH,rcal rad., resummed -^ribn^^ ^2 J ^ ^-L J dk^ 

(- ln(A;, + ze - pK) + ln(/e, + le + (3h)) (22) 
1 2h 



where we have defined the resummation functions, from Ref. [5], 



7, = 2CF^^^^iHs/m^) - 1) (23) 



,2 



2 TT ' 3 2 
and 

Here, Cp is the quark representation quadratic Casimir invariant already defined in (jl]), 
s = {pi + QiY in our process in ([1]) specialized to X = V^*^ with = p\, 

Ce = .5772 . . . 

is Euler's constant and T(w) is Euler's gamma function. The function Fyfs{z) was already 
introduced by Yennie, Frautschi and Suura [11] in their analysis of the IR behavior of 
QED. Using the substitution fc^ = y/ekz, we have 

^7f ^ 2 / / — — 

-^nfen^g—o I lR,real rad., resummed Ffibiv^ j d k± J dk^ 

(- ln(^, + i^e- (3h) + ln(^, + i^e + ph)) (26) 

1 2k_ 

-{I- (3^)ekl-\i^^kl + e' 

We see that the RHS of this last equation vanishes as e 0, removing the violation of 
Bloch-Nordsieck cancellation in ( 120|) . and, thereby, in (I3])|f|. 



^Note that, by the mean value theorem, the RHS of ^ is equal to e'''/'^ Fy Fs{lq)e^''''^ < {Pk^y" > 
-fnbn^ij-o|in,reai rad., IR pole part, whcre < A > dcnotcs thc rcspcctivc mean value of A defined with > 4; 
thus, ((26|) is still a higher twist effect with a coefficient which vanishes as e — > 0. 
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We conclude that the result in Refs. [2] is obviated by amplitude based exact resum- 
mation of the higher order corrections in QCD perturbation theory. Only the infrared 
singular term from (IT^ is exponentiated, so that the finite non-zero terms in the cross 
section are all treated on equal footing - there is then no scheme dependence introduced 
by our resummation. The way is open to employ the current quark masses in ISR phe- 
nomenology for the LHC For the light quarks, their main use will be as collinear/IR 
regulators, as the usual factorization methods [12] will generally replace them with the 
scale of such factorization; for the b quark, we can not exclude at this time that its mass 
may have some additional role in precision LHC theory. Indeed, in addition to current 
algebra constraints, we know that from the measured differences between the parton den- 
sities for s, c, and b quarks in the proton that the "heavy" quark masses can not actually 
be zero. We follow Ref. [130 in defining parton densities for heavy quarks here. The issue 
then is the accuracy of the massless approximation in the ISR in the context of precision 
LHC physics; for, already in QED, it is known that the corresponding limit rrie i in ISR 
and the condition rrie = in ISR differ in 0{a/7r). In QCD ag/vr = 3% at TeV scales and 
this would be unacceptable if it would occur when the precision tag is 1%, as it will be 
at the LHC for some processes. 

We note here that there is considerable literature [14-19] on the use of quark masses in 
perturbative QCD phenomenology, especially for deep inelastic scattering (DIS) processes. 
While in the original ACOT [14] variable flavor number scheme and in Ref. [18], quark 
masses are retained in the initial state analysis, in most cases, following the S-ACOT [17] 
variable flavor number scheme and various extensions [16, 19], the ISR is treated with zero 
quark mass in the hard scattering coefficient with possible use an appropriate rescaling 
variable x(H-4m^/Q^) [15], in standard DIS notation. These anaylses result in general in 
a better fit to the available structure function data, although for Ref. [16] the significance 
of the attendant improved is within the range of uncertainty of the respective fully 
massless result. These efforts all speak to the need for proper treatment of quark mass 
effects in precision high energy QCD phenomenology. 

We have discussed the theorem in Refs. [2] in which the Drell-Yan process for quark- 
quark scattering is considered. However, our solution for the lack of Bloch-Nordsieck 
cancellation only depended on the external lines in the initial state, so it will carry-over to 
all such ISR configurations: exponentiation of real corrections will render an extra factor 
of /cg' in the respective integral over phase-space to remove any end-point contributions 
which are not already canceled by virtual corrections as required by the Bloch-Nordsieck 
theorem. 

Further implications of the results in this paper will appear elsewhere. [20] 

'^In the proof of factorization presented in Ref. [13] for heavy quarks, there is an imphcit use of the 
cancellation of ISR infrared singularities; our results remove any issues concerning this use. 
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Figure Captions 



Figure 1. Graphs evaluated in Ref. [2] (see the first paper therein, especially) in arriving 
at the result inQ using Mueller's theorem for the respective cross section. Here, the usual 
Landau-Bjorken-Cutkosky(LBC) [9] rules obtain so that a slash puts the line on-shell and 
a dash changes the ie prescription; and graphs that have canceled or whose contributions 
are implied by those in the figure are not shown explicitly. 

Figure 2. The contour C used in the complex /c^-plane to evaluate the real emission part 
of the contribution of diagrams (q-o) in Fig. 1 to the RHS of ([3]). See the text for further 
discussion. 
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Figure 1 . Graphs evaluated in Ref. [2] (see the first paper therein especially) in arriving at 
the result in (3) using Mueller's theorem for the respective cross section. The usual 
Landau-Bjorken-Cutkosky (LBC) [10] rules obtain so that a slash puts the line on-shell 
and a dash changes the is-prescription; and, graphs that have cancelled or whose 
contributions are implied by those in the figure are not shown explicitly. 
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Figure 2. The contour C used in the complex kz -plane to evaluate the real emission part 
of the contribution of diagrams (q-o) in Fig. 1 to the RHS of (3). See the text for further 
discussion. 
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